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ABSTRACT 


The paper demonstrates the use of simplified elasticity solutions 
to determine the mechanical response of thick laminated beams and plates 
subject to out-of-plane loading. 

Excellent results were obtained which compare very favorably with 
theoretical numerical and experimental analysis from other sources. 

The most important characteristic of the solution methodology 
presented is that it combines great mathematical precision with 
simplicity. This symbiosis has been sorely needed for design with 
advanced composite materials. 

KEYWORDS: Laminated beam, laminated plate bending, simplified 


elasticity solution 



1. Introduction 


The number of applications for which polymer based composites are 
used is increasing every year. Projections for the last decade of our 
century indicate that structures will be composed of between 25% and 50% 
of composite materials depending on the industrial sector for which they 
are designed and built. The latest developments in the areas of 
processing and manufacturing technology open such a broad avenue of new 
possibilities that the use of composites, by the end of the nineties, 
will probably be even higher than the most optimistic predictions. 

With this increasing number of applications, at least two problems 
become more and more urgent: 

a) How to properly design with composites 

b) How to ensure the durability of components subject 
to mechanical and environmental loading conditions 

The classical laminated plate theories (CPT) are inadequate for 
design especially in those applications where out of plane loads need to 
be supported, i.e., bending, buckling, impact etc. The shear- 
deformation effects are unaccounted for in the CPT which in turn leads 
to incorrect estimates of durability. 

The objective of this paper is to demonstrate the use of a 
simplified theory of elasticity approach for the stress analysis of 
laminated beams and plates. The advantage of this solution methodology 


i s that 

closed 

form 

analytical 

solutions can 

be obtained 

for 

any 

laminate 

layup. 

In 

addition. 

factors which 

are related 

to 

the 


durability such as resin rich regions, gradients in material properties 
due to processing of thick laminates, etc. can be incorporated without 
undue complications. 



The theory of elasticity for laminated beams and plates is based on 
an approach which has been advocated by Biot [1,2]. The accuracy of the 
resulting solutions to determine the deflections, stresses and strains 
in multilayered components is demonstrated. Results are compared with 
existing elasticity solutions, experiments and finite element results. 

The ability to use analytical solutions which are simple, accurate 
and which can incorporate realistic fabrication features should greatly 
enhance the engineers ability to design durable and trouble-free 
laminated components. 


2. Formulation of the Problem 

We consider an elastic and orthotropic plate of thickness h, 
described by the displacement components u and w. The x-axis is 
directed along the span and the z-axis is perpendicular to the plate, as 
shown in Fig. 1. The strain components are: 

_ 3_U _ 3W _1 ,3JJ 3W\ / 1 v 

e x ” ax ’ e z ~ az * e xz ” 2 Sz ax-* ' 

The stress-strain relations for the orthotropic material are: 

a x = C 11 e x + C 13 e z 
a z * C 13 e x + C 33 e z 
T xz = 2C 55 e xz 

in which C u = (1 - v 2 3 E 33 /E 22 ) E n /V 

C 13 * (v 13 + v 12 + "23) E 33 /V 
C 33 = t 1 ‘ E 33 /E ll) E 22 /V 


( 2 ) 


( 3 ) 


2 



C 55 ' G 13 


V = [1 - vp(\i 
2 


12^12 ^22 /E U + 2v 23 v 13 l '33 / "ll 


E -?^/ E ii) " v i3 E 33^11 


- v 23 e 33 /e 22 J 


(4) 


If the stress o z , normal to the plate, is neglected and assumed to be 
zero, then relocations (2) can be rearranged and become. 


\ m 4M e x 
T xz * 2 G 13 e x2 


The coefficient M can be written as: 


M = 4C 33 (C 11 C 33 ' C 13* 


Substituting the relations (3) into (6), we obtain 
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M ' 4 U ' "13 *3? 


and for an isotropic material; 


M = 


4(l-v 2 ) 


(5) 


( 6 ) 


(7) 


( 8 ) 


The equilibrium equations are: 


3 



( 9 ) 


3a 


x 


3 


X 



= 0 


3t 


xz 


3 


X 



- 0 


Note that the assumption o z = 0 is only introduced in the stress-strain 
relations and not in the equilibrium equations. This is due to the fact 
that we can neglect the magnitude of a with respect to a and x but 

Z A a4 

not its spatial rate of change. We investigate displacement fields, 
which are sinusoidally distributed along x. Trial solutions for u and w 
can be written as: 


u = U(z) sin (ax) (10) 
w = W cos ( ix ) 

In which a is the wavelength along the x-di recti on. An additional 
approximation is introduced here by assuming W to be a constant, equal 
to the average displacement across the thickness. The second equation 
(3) yields 


T xz = sin ax 
with x (z) = G 13 - aw) 


By eliminating and U between equations (3) and the first of equations 
(9). We btain the following differential equation in t . 

A 4 


d_ ( 1_ dt xz \ 1 r x 

dz MM dz ' " 


a 3 w 


(12) 


4 



The integration of this second-order differential equation in t xz 
requires two boundary conditions. Let the boundaries of the beam be 
located at z = ± h/2 (Fig. 1) and assume the shear = t xz (h/2) and 
= t (- h/2) to be given at the top and bottom. With these boundary 
conditions the function t (z) is obtained by integration of Eq. 12, 

A 4 

where M, G^, i and w play the role of parameters. By integrating the 
second equilibrium equation (9) along z, we obtain: 

h/2 

q - - J t x _ dz (13) 

-h/2 xz 

The total load applied to the same unit area is: 


( a 2 1 i = [o z \ 2 = 9 COS1X 


(14) 


Since t(z) is known in terms of w, equation (13) determines the 
deflection w when the load q is given when we know t , the values of U 

A4 

and o x are determined by combining the first of equations (1) with the 
first of equations (3). We obtain: 


U = 


xz 


4M& 


2 dz 


_ J. dT xz „ eBV 
°x ■ i dz C0S!,X 


(15) 

(16) 
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3. Solution for Lamina and Laminates 
The foregoing methodology may be used to analyze laminates 
constituted by the superposition of adherent homogeneous layers. 

Consider first a single-orthotropic lamina of thickness h and 
constant elastic coefficients and M. The shear stresses at the top 
and bottom of the lamina are denoted by x^ and t 2 respectively. 

Eqn. 12 is readily integrated in this case, the solution of which 
is 


t xz = cosheaz + C 2 sinhsaz - i G 13 W 


(17) 


with 


c i " ( T i + t 2^ + 1 


13 1 cosher 


r _ 1 , s 1 

~ 2 lT z ' t 2' sin 


nhBy 


1 

2 


y = — ih 


B » 2/ 


’13 


From equation (15) we derive the values and U 2 of U at the top and 
bottom of the layer 


U l = 


1 


4 / MG 


13 


U 2 = 


1 


4 / MG 


(x^a + x 2 b) + cW 


(x^b + x 2 a) - cW 


( 18 ) 


13 


with 


a = tanhBY + 


1 


tanhBY 


6 


b ’ tanhs ’' - tiFKiT 


c = -^ tanhSY 
0 

The normal load q cosax applied to this layer is obtained from Eq. 13. 

q - - (tj + t 2 ) c ♦ A g 13 w(i - £) (19) 

For a laminate constituted of N adhering layers, the i-th layer of 
thickness h^ is characterized by the material properties 

(Gj 3 ) and (M 1 ) . Corresponding parameters are a^, b^ and c^. We denote 
by t and t . + ^ the shear stresses at the top and bottom of the i-th 
layer respectively and by and the displacements at the 

corresponding faces. The condition of adherence of the layers i and i+1 
are obtained by equating the displacements at the interface. Applying 
equations (18), we obtain a recurrency equation which has to be 

satisfied by the interface stress at three subsequent interfaces. 


B, ,, 


(A, 


+ * 1+1 


> r i+l 


’ 1+1 


r i=2 ■ ' (C i 


C 1+l> 


iW 


i 

where A. * — — 1 — i = 1,2 ... N 

4/ M 1 g{ 3 


4/ M 1 G{ 3 

■1 _ G LT G TT 

1 T ~ ? ? 

G, T sin e + G tt cos e 

Ll 1 1 layer # i 


( 20 ) 


(21) 
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layer # i 


( 22 ) 


M 1 


’ll 


4 0 - f^“) 

u t u 


with: 


’ll 


'13 


1 

COS 2 0 . , 1 2vLT\ ..2 „„_2, 

r (p " ri ) sin 9 COS 6 

e l V g lt el 

E U (^I (sin 4 0 + cos 4 e) - 


sin 2 e 


111 2 ? 
(•?- + -F — «-=-) sin 0 cos 0 l 

fc L t T b LT 


note: for 0=0 


Thus 


for 0 = 90° 


Thus 


G 13 = G LT = G 13 


E 11 " E L " E 11 


v 13 3 y LT * y 13 


M = 


’ll 


2 E 2, 


4 d - v 13 


G 13 = G TT = G 23 


E 11 = E T = E 22 


E T E 

v 13 ' \T E, 3 v 13 E 


M = 


•22 


4 d - vf 3 


22 

11 
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in which 9 is the angle between the fibers (which are oriented in the x- 
y plane) and the x-direction. 


Eqn. 20 is a system of (n-1) equations, in the (n-1) interface 
stresses. Assuming that the interface stresses and t n at the top and 
bottom of the laminate are given. The other interface stresses are 
obtained in terms of a single unknown W. The latter is evaluated by 
considering the total load of q cos ix applied to the laminate. It is 
the sum of all individual loads of cos ax acting on each layer. Hence 

q - l q i (23) 


where according to Eq. 19 


C, 


"i - (V *i*i> S * 1 h i G « < 1 - 77 > w 


we may write 


Vl> c,- + A I h i sis' 1 - 77 


1 

c i 

1 


(24) 

(25) 


Since t. is a known function of W, while q is given, equation (25) 
determines W. 

Equations (13 ... 21) thus allow us to obtain closed form 

analytical solutions for any laminate layup. These solutions can easily 
be programmed on a hand calculator or PC and be used in the design 
office for preliminary sizing of laminated composite beams and sandwich 
panels. 
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In the remainder of this paper, we will demonstrate the use and the 
accuracy of this approach. Results are compared with theoretical, 
numerical and experimental data which is available in the literature. 


4. Results for Lamina with 0° and 90° Orientation 

We shall first illustrate that well-known strength of materials 
results can be retrieved or a special case of the elastcity solution 
which we outlined in 3. 

Consider the case of a simply supported homogeneous orthotropic 
beam subject to three point bending. The conditions at the top and 
bottom surface are = 0 and = 0, as schematically indicated in Fig. 
2 . 

Upon substitution of and in Eg. 19, we obtain 


q - A E 13 (1 - ****») w 


(26) 


or W 


a 

a 2 h (1 


tanhsy ) 

3y 


(27) 


The deflection of a beam on three-point bending is easily found by 
expanding the load P in Fourier series and by adding the corresponding 
displacement components. The result is: 


,, _ 2P y n 1 

aa 2 hG n 1,3,5 (1 - tanhs-^T 

% 


(28) 


with Y m = m 


Th 

2a 
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Expanding tanh8Y m in series, and retaining the first three terms, we 
obtain 


PE 


W = 


11 


2 * < S 13> 2 


1 ” 5 |< 8 I 1 'IS <%> 2 ' 


(29) 


using the binomial approximation 


1 


6 

15 


1 - (sr m )' 


6 

15 


* 1 + * O’.)' 


(30) 


we obtain, by adding the three terms 


W = 


Pa' 


48 E 11 1 


E n h 2 
[.9994 + 1.12 — (?) 1 

8 13 


(31) 


This is a well-known result, from the theory of strength of materials, 
that is often used to obtain the interlaminar shear modulus (G^) from 
short beam deflection data. 

The difference (in %) between the deflection obtained using the 
elasticity solution (Eq. 28) and the strength of materials result Eq. 
31) is plotted in Fig. 3. The curves are representative for a 


undirectional graphite epoxy with E^ = E^ * 171 Gp a (25 10° psi) 
8 13 = 8 13 3 Gpa (*210 psi) ♦ The difference is calculated as: 


and 


DIFF = M (Ed- 31) - W (Eg. 28) im 
uirr W (Eq. 28) x iUU7 ° 


(32) 


Fig. 3 indicates that the strength of materials solution slightly 
overestimates the deflection (less than 2%) for ratio's of beam-span to 
thickness which are larger than 10. Both solutions deviate rapidly for 
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span to thickness ratio's which are less than 4. We can thus conclude 
that the strength of materials formula (Eq. 21) is a good approximation 
for all practical purposes. 

Unfortunately simple equations as Eq. 31 are not available for 
laminates: Nevertheless in 5 we will present actual solutions obtained, 
using the methodology outlined in 4. These solutions represent an 
excellent compromise between precision and flexibility as will be shown. 


5. Results for Laminates 
5.1 Cylindrical Bending of a 0°/90 o /0° Laminate 

The solution methodology for a laminate is first applied to a 
0°/90°/0° subject to cylindrical bending, as schematically shown in Fig. 
4. The results obtained are compared with these which were published by 
Pagano [3] who rigorously solved this elasticity problem (i.e. without 
bringing in the assumption that a z is much smaller than a x and t ). 

First Eq.'s 20 are solved for the shear stresses at the 0°/90° 
interfaces. Upon substitution of these results into Eq. 25, the center 
deflection is obtained as 


W = 


i [- 


/ Mj_ G 13 


2 b 2 C 2 
/ m 2 g 23 


+ 2ih^ (1 


''I 

-7) + ih 
Y 1 


2 G 23 


(1 -7^ 

y 2 


where coefficients a^, b^ and c^ are defined in Eq. 18 and and M 2 are 
defined by substituting 0 = 0° and 90° into Eq. 22. The same 
substitution into Eq. 21 defines G^ = G L y and G 2 3 = Gjj. 
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The numerical values used Pagano's elasticity solution are: 


E n = 171 GPa (25 10 6 psi) 
E 22 = 6.85 GPa (10 6 psi) 
G 13 = 3.42 GPa (.510 6 psi) 
G 23 = 1.37 GPa ( . 210 6 psi) 


The plate has unit thickness, consequently h i = h 2 = I ( as ideated in 
Fig. 4). 

In order to compare our results with Pagano's [3]. We need to 
multiply the deflections, obtained from Eq. 33, with a normalization 
factor: 


W = W 


100 E 22 (h 2 + 2h 1 ) 


(34) 


V 


The value W is obtained for different span to thickness ratio's (a/h) 
and represented by the solid line in Fig. 5. The circles are 
the W values which were obtained, using the much more cumbersome Pagano 
solution. As can be seen, the agreement is excellent. Eq. 33 is thus 
especially useful as a simple solution for short thick beams (as used 
for composite leaf springs, hinges, helicopter rotor-connections and 
many other primary loadbearing structural applications). Fig 5 also 
indicates that our solution for W correctly converges towards the 
classical plate theory (CPT) solutions for large span to thickness 
ratios. 

The shear stress t 2 was solved, from Eq's 20, as: 


t 2 


= AW a 


(35) 


13 



with A = - 


(35) 


4(C 1 + C 2 ) 


/ G 13 


2 ^2 

/ M 2 g 23 


To obtain this solution for t 2 we made use of the symmetry of the 
laminate, with respect to the x-axis, thus t 2 = t 3 and = 0 . 

After substitution of Eq. 35 into Eq. 17, we obtain the shear 
stress distribution in the 0° layer (layer 1) 


T xz = C 1 cosh8 l 12 + C 2 sinh8 2 12 ~ 1 G 13 W 


with C. = [iy A + G, 0 ] 


2.W 


13 J cosh 8^ y ^ 


(36) 


r _ 1_ A tw 


'2 2 sinh 8j_ y^ 


The shear stress distribution in the 90° layer (layer 2) is 
obtained as 


t xz = cosh8 2 iz - iG 23 W 

(37) 

with C. = (A + G«~) 

1 v 23' cosh b 3 Y 2 

The distribution of the interlaminar shear stress through the 

thickness can now be plotted and compared to results obtained by Pagano 
[3], as shown in Fig. 6. The solid line is the current solution, as 

represented by Eq. 1 s 36 and 37 respectively, whereas the circles 

represent the solution obtained by Pagano. Both solutions were obtained 
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for a span to thickness ratio of 4. As can be seen, the agreement 
between both solutions is excellent. The shear stress reaches a maximum 
at points A and C (Fig. 6.) in the 0° plies, which leads to the well- 
known cusp-like features in the shear stress distribution. The 
mathematical criterion for the occurence of a maximum within the 0° ply 
can be obtained by putting the derivative of Eq. 36 with respect to z 
equal to zero and by solving this equation for z. This procedure is 
demonstrated in appendix A. 

It is very encouraging that the current solution-methodology 
enables us to obtain the detailed features of the interlaminar shear 
stress distribution. The interlaminar shear stress is matrix dominated 
and has a very damaging effect on the integrity of structural composite 
components. The obtained result is essential since this work is part of 
a larger program on durability of composites and wishes to incorporate 
the time-dependent aspects which are due to the matrix-dominated nature 
of these stresses. 

Additionally, it has been shown that when tough thermoplastic-resin 
composites are impacted, most of the energy is dissipated through 
interlaminar shear deformtaion [ 4 J . This again emphasizes the 
importance of this solution methodology for the design of impact 
resistant structures. 

The distribution of bending stresses in the 0°/90°/0° laminate is 
obtained, using Eq. 16 combined with Eq. 36 (for the 0° ply) and Eq. 37 
(for the 90° ply). The solid line in Fig. 7 was obtained using both 
formulas, whereas the circles represent the Pagano solution [3]. 
Excellent agreement is obtained. 
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By comparing Fig's 6 and 7 it can be seen that the corresponding 
locations, labeled A, B, C, in both figures coincide with two maxima and 
minima in the shear stress distribution. It is well known from the 
theory of strength of materials, that the bending stress at the neutral 
axis (point B) is zero, while the shear stress reaches a maximum. 
Strength of materials does not predict more than one maximum and is thus 
inapplicable for complex-multilayered composites which might have a 
number of maxima, proportional to the number of plies. 

5.2 Modeling of Interleaves or resin-Rich Regions 

Current methods to reduce the susceptibility of laminates to impact 
damage rely on the use of adhesive layers as interleaves . It was found 
[51 that these suppress impact- induced delamination, toughen the 
interface between two lamina and reduce matrix cracking. 

The methodology outlines in this paper should serve well to guide 
in the optimal placement of adhesive layers to obtain the greatest 
benefit in improving impact resistance. 

Consider a laminate, which consists of four 0° - plies, being 
subject to cylindrical bending. The loading conditions, geometry and 
material properties are same as for the example given in 5.1. 

The interlaminar shear stress-distribution in the four-ply laminate 
was calculated using Eq. 17. The result is represented by the solid 
line in Fig. 8, whereas the circles represent results obtained by using 
Pagano's [3] solution methodology. 


16 



The introduction of three epoxy interleaves (E = .68 GPa or 10^ psi 
and \> = .3) each of which only represents 3% of the total thickness 
drastically alters the shear-stress distribution as represented by the 
dotted line in Fig. 8. 

As can be seen, the maximum shear stress, reached at the center of 
the laminate, is 33% higher than for the case of the ideal 4-ply 
composite. This might exhaust the shear-deformation capability of the 
epoxy and lead to failure. Interleaves are thus more effective to 
suppress impact when placed closer to the surface. The presence of 
resin rich area's also enhances the creep-characteristics of the 
laminate and reduces its buckling resistance. Both these features are 
currently being implemented into the methodology which we outlined in 2. 

It can also be seen in Fig. 8 that the shear stress remains 
constant through the thickness of the interleaves. Consequently there 
will be discontinuities in the shear strains, because for shear modulus 
of the adhesive we used 1.37 GPa (.2 lb 6 psi), while the shear modulus 
of the plies was taken as 3.42 GPa (.510 6 psi). 

5.3 Experimental-Theoretical Correlation for Thick Laminates Subject to 

Three-Point Bending 

Post and coworkers [6] obtained interlaminar shear strain 
distributions on 48 ply laminates subject to three-point bending. This 
experimental data was obtained by means of a high-precision mo ire 1 - 
interferometry method. A schematic of the experimental setup and the 
dimensions is schematically shown in Fig. 9. A shear strain 
distribution was experimentally obtained along the line A-B, which is 
located at quarter span. 
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The analytical solution, which we obtained for the shear strain. 


using Eq. 17 can be written as: 


_ 2P 
r xz ~ azh G 


13 


1,3,5 


,cosh_m8£z 

'cosh 8 y_ " ' 

m 


(i - 


tanh 8y 


nu 


sin Y m x 


(38) 


By, 


m 


All variables in Eq. 38 have been defined in section 4. The applied 
load is represented by P * 6120 N (1350 lbs). The thickness h = 12.78 
mm (.5") and the span is 63.5 mm (2.5"). 

The material used was T300/5208, which has the following material 
properties: 

E u = 130 GPa (19 10 6 psi) 

G 13 = 6.85 (10 6 psi) 

The series in Eq. 38 was evaluated for up to 200 terms. A 
comparison revealed that changes of less than 1 % occurred after the 
first eight to ten terms. The resulting shear strain distribution is 
represented by the solid line in Fig. 10. The circles represent the 

experimentally measured strains. As can be seen, the obtained agreement 
is very good. 

This agreement suggests a possible approach to obtain the 
interlaminar shear moduli G^ and G23. First the analytical formula, 
Eq. 38 is used to determine the maximum shear strain at quarter span. 
This maximum shear strain is subsequently plotted for shear moduli 
ranging from 1.37 GPa (.2 MSI) to 13.7 GPa (2 MSI). The result is 
plotted in Fig. 11. It is now sufficient to mark the experimentally 
measured maximum strain on the ordinate and to read the corresponding 
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interlaminar shear modulus, on the abcissa, as indicated by the 
arrows in Fig. 11. The value thus obtained for G^ is 6.23 GPa (.91 
MSI). We should add that this is exactly equal to the value reported 
for T300/5209 by Sumsion and Rajapaske [7], Data obtained by means of 
ultrasound measurements on T300/5208 [ 8 ] is the only other literature 
source we found and which lists a value of 7.12 GPa (1.04 msi) for Gjj. 

The approach should also work to determine the shear modulus G^, 
which is often considered to be equal to G^ by the engineers who 
perform the finite element analysis. If indeed we would have made this 
assumption, to solve the problem for the 0°/9Q°/0° laminate in 
cylindrical bending, (see 5.11) the cusp-like features in the shear 
stress-distribu-tion (Fig. 6 ) would not have shown up, as discussed in 
appendix A. 

The experimental results for T300/5208 and T300/5209 indicate that 
G 23 is 44% and 49% lower than G^. A difference of this magnitude is 
important enough to be considered for various design applications. 
Unfortunately the data is sorely lacking. Some data is available at 
room temperture, while none is available at elevated temperature. This 
conclusion is also true for the case of metal-matrix composites. 

4.4 Load-Deflection Behavior of a Shear-Beam 

Moussiaux et al. [91 recently published a paper on a new short-beam 
shear test specimen. The objective is to use the specimen for the shear 
characterization of an adhesive. The test-setup is schematically shown 
in Fig. 12. 

Two aluminum adherends are bonded together with an adhesive. The 
methodology developed in [91 relates the overall deflection of the beam 
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methodology which we outlined in this paper can be effectively used to 
obtain results which are identical to these in ref [9J. We should add 
that the current approach is more general, because anisotropic adherends 
can be accounted for. It is also possible to locate the adhesive layer 
off-center and thus subject it to combined tensile and shear stresses. 

The deflection at the free end of the beam is obtained by expanding 
the load P in Fourier series and by applying the summation procedure, 
discussed in section 4, to Eq. 33. This equation is in principle only 
applicable for the case of a single supported beam. The current results 
for a clamped beam were obtained by using the statically equivalent case 
of a single supported beam of twice the length and subject to twice the 
load. 

Deflections were calculated for aluminum adherends (with E = 68. 
GPa or lO^psi and v * .3) and for varying properties of the adhesive. 
The length of the beam a = 63.5 mm (2.5"), the thickness of the 
adherends h^ = 6.35 mm (.25") and the thickness of the adhesive was 7.27 
mm (0.5"). The applied load was 445 N (100 lb). 

First we calculated a reference deflection, W ref assuming that the 
adhesive has the same properties as the adherends. (i.e. the deflection 
of a clamped aluminum beam which has a thickness of 13.97mm, (.55")). 

Subsequently the adhesive was made more and more compliant and the 
calculated deflection of the beam consequently increased. this 
deflection was normalized with respect to the reference deflection. 

The ratio W/W re f has been plotted in Fig. 13 in which the abcissa 
is the ratio of young's modulus of the adherends (E) and the shear 
modulus of the adhesive (G n ). The solid line has been obtained using 
the methodology outlined in this paper. It was assumed that the 
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adhesive and the adherends have the same poissons ratio (v = .3). The 
ratio E/Ga thus increases from 2.6, for the solid aluminum beam, to 10^ 
for adherends which are basically no longer bonded. The triangles and 
the square symbols in Fig. 13 were obtained by Moussiaux et al. [9] and 
by using finite elements. 

The deflection treshold, which is obtained for high E/Ga ratio's is 
reached because the very low stiffness of the adhesive allows the two 
adherends to act independently. A beam with unbonded adherends thus 
deflects about five times more than a solid aluminum beam. 

This value for the deflection treshold can also be obtained on the 
basis of strength of materials. The reference deflection is given as: 


W 


ref 


12Pa 3 

3E(2h 1 + h 2 ) 3 


(39) 


The debonded beam can be pictured as two separate beams with a thickness 
h^, each of which is acted upon by a load of P/2. We thus obtain: 


W = 


P a 3 
2 a 


3E 


12 


Dividing (40) by (39) we obtain: 


JL_ s i ( 2 + _i) 

w ref £ n l 


(40) 


(41) 


Upon substitution of the numerical values for h^ and f^, we find 


W 


ref 


= + - 


5.32 


(42) 
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which is very close to the result obtained in Fig. 13. 

Assuming that the adhesive was an epoxy with E = .68 GPa (10^ psi) 
and \> = .3 we obtained the maximum shear stress in the adhesive layer as 
a function of the relative position along the beam, as plotted in Fig. 
14. The circles represent the results obtained by Moussiaux et al. 
[9]. Again the agreement is excellent. 

It is generally agreed upon that the interphase between the 
adherent and the adhesive needs modeling. The solutions which we 
presented are ideally suited to accomplish this task. Interphase 
regions which are only a few hundred angstrom thick can be easily 
accounted for because of our ability to derive analytical solutions 
there is no risk to run into mathematical problems or problems with 
computer algorithms. Precise solutions are guaranteed. 

Results obtained by FIOR and Brinson [ 30 ] indicate that the use of 
deflection measurements on the shear beam cannot be used to determine 
the shear modulus of stiff adhesives. Instead, shear deformation 
measurements similar to the procedure discussed in section 5.3 have been 
recommended. 

Notwhitstanding these comments, we further developed this procedure 
into dynamic regime, such that a simple frequency measurement defines 
the shear modulus of rigid as well as compliant adhesvies. These 
procedures are outlined in [11]. 


5. Conclusions 

We demonstrated the use of simplified elasticity solutions, which 
lead to precise results for the displacements and the stresses in 
laminated composite beams and plates. 
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These solutions can be programmed on a simple hand calculator or on 
a PC. They can be used in the design office for preliminary sizing of 
structural composite and sandwich components subject to out-of-plane 
loads. 

It was shown that excellent agreement was obtained with the much 
more cumbersome full-elasticity solutions for laminated plates. Special 
features, such as adhesive interlayers and resin-rich regions can be 
accounted for without additional complications. Excellent agreement was 
obtained with experimentally measured shear strains, obtained on a 48- 
ply laminated beam which was subject to three-point bending. This 
agreement led to the proposition of a new interlaminar shear modulus 
characterization test. 

It was also shown that good agreement was obtained with results for 
the adhesive bonded shear beam. This opens the way to in-si tu 
measurements of adhesive bond line properties. 

The solutions which we discussed are useful to support and to 
advance a framework for durability predictions that has been proposed by 
Brinson [12]. More realistic constitutive equations, including 
hygrothermomechanical effects can be included by using the internal 
variable approach. 


£ x’ e z ,e xz 

c ll’ c 13’ c 33’ c 55 


W T xz 

u,w 

M 


6. Nomenclature 
strain-tensor components 
coefficients of the stiffness matrix 
stress-tensor components 
displacements in the x- and z-directions 
stress-strain law coefficient (defined in Eq. (6)) 
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EijEj, v 

g lt* g tt 


LT* 


^11 E 33’ v 13* 

b 13’ G 23 


wavelength 

transverse load distribution 

dimensionless anisotropy factor (defined in eq. 17) 

dimensionless geometry factor (defined in eq. 17) 

displacements along the top and bottom of a ply 

interlaminar shear stresses at the top and bottom of 
a ply 

transverse load 
span 

thickness of a ply (or laminate) 
material properties using the L-T notatin 

The same material properties, using the 1-3 notation 
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8 . Appendix A: Maximum Shear Stress in a 0°/90°/0° Laminate 

As can be seen in Fig. 6 , the shear stress reaches a peak in the 0° 
layer. Eq. 36 can now be used to derive an expression for the location 
of this peak the shear stress is maximum if 

H-0 (A. 1) 

32 V 

Thus after subsitution of Eq. 36 into (A.l) we obtain 


sinh 8 ^ %z + C z cosh iz * 0 


(A. 2) 


Solution of A. 2 for z and substitution of the values for and C 2 
defined in eq. 36 gives 


z = tanh - ^ [ 

8 ^ 2 . 


1 + 


2i G 
t 2 


13 


•] (A. 3) 


W) tanh 


hi h^ 

Note that z lies in the domain — ^ to + j- . 

The obtained value of z is a local coordinate measured to reference 
axes which are located at the center of the 0° ply, as indicated in Fig. 
A.l. 

Eq. A. 3 was used to determine the maximum shear stress in the 
0°/90°/0° laminate, for a varying G 23 modulus of the 90° layer. All 
other data used was the same as given in example 5.1. It was assumed 
that the laminate has unit thickness (thus each ply has a thickness of 
1/3). Table A.l contains the obtained results. 
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The negative signs in the solution for z mean that the maximum 

shear stress develops in that part of the 0° plies closer to the 0°/90° 

interface. It can be seen in the table that an increase in shear 

modulus G23 drives the maximum into the direction of the interface. 

This is made clearer in Fig. A. 2, where we normalized G 2 3 with respect 

to G|3 (G^ = 5.10® psi). The vertical axis identifies the location 

h l 

of (t ) max in % of — ~ (100% means the maximum is at the 0°/90° 

aZ l 

interface, 0% means the maximum is located at the center of the 0° ply). 

As can be seen from the figure, the maximum shear stress will act 
at the interface for G23/G ^3 ratio's larger than 7. A very low G23/G13 
ratio will force the 0° plies to act independently and force the maximum 
shear action at the center of each ply. The same reasoning applied to 
adhesively bonded beams, in which the adhesive changes properties due to 
environmental degradation and/or viscoelastic effects. 
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List of Captions 


1. Orthotropic composite plate 

2. Simple-supported orthotropic beam 

3. Deviation between deflections obtained with strength of 
materials formulas and formulas outlined in this paper (for 
graphite epoxy material) 

4. Cylindrical bending of 0°/90°/0° laminate. 

5. Comparison of deflections for a range of span/thickness ratio's 
obtained using the current model (solid line) and Pagano's full 
elasticity solution (circles). 

6. Comparison of interlaminar shear stresses, for span/thickness 
4, obtained using the current model (solid line) and Pagano's 
full elasticity solution (circles). 

7. Comparison of bending strains, for span/thickness = 4, obtained 
using the current model (solid line) and Pagano's full 
elasticity solution (circles). 

8. Effect of interleaves on the interlaminar shear stress. Curve 
1 line = no interleaves. Curve 2 = three interleaves. 

9. Experimental setup and dimensions of a 48 ply 0° - deg 
T300/5208 laminate 

10. Comparison of experimental (circles) and analytically (solid 
line) obtained shear strain distributions at quarter span. 

11. Variation of maximum shear strain, at quarter span, as a 
function of the shear modulus. 

12. Setup and dimensions of short beam shear test specimen. 

13. Deflection amplification for bonded shear beam for increasingly 
compliant adhesive. 

14. Maximum shear stress in the bonded shear beam for various 
positions along the span. Comparison with solution by 
Moussiaux et al. (circles) 


28 



Appendix B 


A.l Detail of 0°-ply local coordinates 

A. 2 Position of the location of maximum interlaminar shear stress 
in the 0° ply 

Table A.l. 
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